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Abstract 



We present a new, fully covariant and manifestly gauge-invariant expres- 
sion for the temperature anisotropy in the cosmic microwave background ra- 
diation resulting from scalar perturbations. We pay particular attention to 
gauge issues such as the definition of the temperature perturbation and the 
placing of the last scattering surface. In the instantaneous recombination ap- 
proximation, the expression may be integrated up to a Rees-Sciama term for 
\/~} [ arbitrary matter descriptions in flat, open and closed universes. We discuss 

the interpretation of our result in the baryon-dominated limit using numer- 
ical solutions for conditions on the last scattering surface, and confirm that 
for adiabatic perturbations the dominant contribution to the anisotropy on 
intermediate scales (the location of the Doppler peaks) may be understood in 
& • terms of the spatial inhomogeneity of the radiation temperature in the baryon 

rest frame. Finally, we show how this term enters the usual Sachs- Wolfe type 



calculations (it is rarely seen in such analyses) when subtle gauge effects at 
■ the last scattering surface are treated correctly. 

1 Introduction 

The calculation of the primary temperature anisotropy in the cosmic microwave 
background radiation (CMB) resulting from density perturbations has a long his- 
tory, beginning with the seminal paper by Sachs and Wolfe M. Since the original 
Sachs- Wolfe estimate, a wealth of detailed predictions for the anisotropies expected 
in various cosmological models have been worked out. The calculations are straight- 
forward in principle, but, like many topics in cosmological perturbation theory, are 
plagued by subtle gauge issues 0. 

The problems of gauge-mode solutions to the linear perturbation equations and 
the gauge-ambiguity of initial conditions can be eliminated by working exclusively 
with gauge-invariant variables, as in the widely used Bardeen approach and the 
less well known covariant approach advocated by Ellis and coworkers ]§, g]. However, 
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gauge issues still arise in connection with the definition of the temperature pertur- 
bation and the placement of the last scattering surface 0, ||. The latter gauge issues 
do not arise at first-order in numerical calculations which integrate the Boltzmann 
equation in a perturbed universe, since the visibility function (which determines the 
position of the last scattering surface) multiplies first-order variables giving only a 
second-order error from the use of a zero-order approximation to the visibility 0. 
However, this is not always the case in Sachs- Wolfe style analyses, which integrate 
along null geodesies back to the surface of last scattering, unless care is taken to 
ensure that the final result involves only first-order variables on the last scattering 
surface, which then only need be located to zero-order. 

In this paper, we present a new expression for the CMB temperature anisotropy 
arising from linear scalar perturbations which is fully covariant and manifestly gauge- 
invariant. We obtain our expression by integrating the covariant and gauge- invariant 
Boltzmann equation []7|, [5] along observational null geodesies, paying careful attention 
to the gauge issues discussed above. Unlike some covariant results in the literature 
(see, for example, @, £§), the expression derived here can be integrated trivially, in 
the instantaneous recombination approximation, up to a Rees-Sciama term in uni- 
verses with arbitrary matter descriptions. (The covariant results in || |9|] can only 
be integrated in baryon-dominated universes, thus excluding CDM dominated uni- 
verses, and other such models favoured by observation.) We base our treatment on 
the physically appealing covariant and gauge-invariant formulation of perturbation 
theory, as described in |4], |5|. In this approach, one works exclusively with gauge- 
invariant variables which are covariantly-defined and hence physically observable in 
principle. The covariant method has many advantages over other gauge-invariant 
approaches (such as that formulated by Bardeen ||). Most notably, the covariant 
variables have transparent physical definitions which ensures that predictions are 
always straightforward to interpret physically. Other advantages include the unified 
treatment of scalar, vector and tensor modes, a systematic linearisation procedure 
which can be extended to consider higher-order effects (the covariant variables are 
exactly gauge-invariant, independent of any perturbative expansion), and the ability 
to linearise about a variety of background models, such as Friedmann- Robert son- 
Walker (FRW) or Bianchi models. 

For universes which are baryon-dominated at last scattering, our expression for 
the temperature anisotropy may be compared to other gauge-invariant analytic re- 
sults in the literature. We show that, with suitable approximations, the result de- 
rived here reduces to that given by Panek |10| and corrects a similar result given by 
Dunsby recently 0. For the baryon-dominated universe, we use numerical results for 
the covariant, gauge-invariant variables on the last scattering surface, obtained from 
a gauge- invariant Boltzmann code , to discuss the different physical contributions 
to the primary temperature anisotropy. In particular, we show that on intermediate 
and small scales, the "monopole" contribution to the temperature anisotropy is de- 
scribed by the spatial gradient of the photon energy density, in the energy-frame, on 
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the last scattering surface. Since the (real) last scattering surface is approximately a 
surface of constant radiation temperature (so that recombination does occur there), 
the inhomogeneity of the radiation energy density in the energy-frame determines 
a distortion of the last scattering surface relative to the surfaces of simultaneity in 
the energy-frame. The extra redshift (due to the expansion of the universe) which 
the photons incur due to the distortion is seen as a "monopole" contribution to 
the temperature anisotropy on intermediate scales. There is a significant "dipole" 
contribution to the anisotropy on intermediate and small scales, which we discuss 
also. 

We end with a discussion of the gauge issues inherent in the original Sachs- Wolfe 
calculation of the CMB anisotropy HJ, focusing on the "monopole" contribution to 
the temperature anisotropy on intermediate scales, described above. This contribu- 
tion is often missed in Sachs- Wolfe type calculations through an incorrect treatment 
of gauge effects at the last scattering surface |J. (Equivalently, the term is of- 
ten missed through a failure to recognise the direction-dependence of the "expected 
temperature" used to define the temperature perturbation in many calculations.) 
This often neglected term, which is not important on large scales, is an essential 
component of the Doppler peaks in the CMB power spectrum. 

We employ standard general relativity and use a (H ) metric signature. Our 

conventions for the Riemann and Ricci tensors are fixed by [V a , Vh]w c = —TZ a bd C u d , 
and IZab = 1Z a bc . We use units with c = G — 1 throughout. 

2 Covariant Cosmological Perturbations 

In this section, we summarise the covariant approach to perturbations in cosmol- 
ogy B HI to establish our notation and conventions. We begin by choosing a velocity 
u a , which is defined physically in such a manner that if the universe is exactly FRW 
the velocity reduces to that of the fundamental observers. This property of u a is 
necessary to ensure gauge-invariance of the variables defined below. We refer to the 
choice of velocity as a frame choice. For most of this paper, it will not be necessary 
to make a frame choice. From the velocity u a , we construct a projection tensor 
into the space perpendicular to u a (the instantaneous rest space of observers whose 
velocity is u a )\ 



where g ab is the metric of spacetime. We use the symmetric tensor h ab to define a 
spatial covariant derivative ( 3 ^V a which acting on a tensor T h '" c d...e returns a tensor 
which is orthogonal to u a on every index: 



Kb = 9ab ~ u a u b , 



(2.1) 




d...e 



= h a p h b r ...hth t d ...h^T r - s 



t...U1 



(2.2) 



where V a denotes the usual covariant derivative. 



3 



The covariant derivative of the velocity decomposes as 

V a Ub = -coab + cr ab + \6h ab + u a w b , (2.3) 

where zu a b = ^[ab] is the vorticity, which satisfies u a zu a b = 0, <7 a b = is the shear, 
which is orthogonal to u a and traceless, = V a u a = 3H measures the volume expan- 
sion rate (H is the local Hubble parameter), and w a = u b V b u a is the acceleration. 
In an exact FRW universe the vorticity, shear and acceleration vanish identically. 
We regard them as first-order variables (denoted 0(1)) in an almost FRW universe, 
so that products of such variables may be dropped in expressions in the linearised 
calculations we consider here. Other first-order variables may be obtained by tak- 
ing the spatial gradient of scalar quantities. Such quantities are gauge-invariant by 
construction since they vanish identically in an exact FRW universe. We shall make 
use of the co moving fractional spatial gradient of the density of a species i, 

X® = ^ (3) V aP «, (2.4) 

and the comoving spatial gradient of the expansion 

Z a = S^V a 9. (2.5) 

The scalar S is a local scale factor satisfying 

S = u a V a S = HS, {3) V a S = 0{l), (2.6) 

which removes the effects of the expansion from the spatial gradients defined above. 
The vector X® is a manifestly covariant and gauge-invariant characterisation of the 
density inhomogeneity. 

The matter stress-energy tensor T ab decomposes with respect to u a as 

T a b = pU a U b + 2u {a q b) - ph ab + 7T ab , (2.7) 

where p = T a bU a u b is the density of matter (measured by a comoving observer), 
q a ee h h a T bc u c is the energy (or heat) flux and is orthogonal to u a , p = — h ab T ab /3 is 
the isotropic pressure, and the symmetric traceless tensor 7i ab = h c a hfT C( i + ph ab is 
the anisotropic stress, which is also orthogonal to u a . In an exact FRW universe, 
isotropy restricts T ah to perfect-fluid form, so that in an almost FRW universe the 
heat flux and isotropic stress may be treated as first-order variables. 

The photons are described by a covariant distribution function f^(E, e), where 
E = p a u a is the energy of a photon with momentum p a , and e a is unit spacelike 
vector along the direction of propagation in the frame defined by u a . The photon 
energy density p^\ the heat flux and the anisotropic stress ir^ are given by 
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integrals of the three lowest angular moments of the distribution function: 



,(7) 



71 



(7) 
ab 



dEdQE 3 f^{E,e) (2.8) 
J dEdQE 3 f^{E,e)e a (2.9) 
J dEdnE 3 f^(E je )e a e b + \p^h ab , (2.10) 



where dfl denotes an integration over solid angles. Higher-order symmetric traceless 
spatial tensors can be used to characterise higher moments of the distribution func- 
tion (see, for example, [|ll]]), and are useful in numerical simulations of the CMB 
anisotropy |7j. We use the temperature difference from the mean (the full sky aver- 
age) as our definition of the temperature anisotropy &r(e), so that 



A8 T {e) = ^f dEE 3 f^(E,e) - 1. (2.11) 



The temperature perturbation 5x{e) is covariantly defined and gauge- invariant (it 
vanishes in an exact FRW universe), and is observable directly. This should be con- 
trasted with the gauge-dependent temperature perturbation used by some authors 
(see Section |5] for examples) . 

The final first-order gauge-invariant variables we require derive from the Weyl 
tensor W a b c d, which vanishes in an exact FRW universe due to isotropy and homo- 
geneity. The electric and magnetic parts of the Weyl tensor, denoted by 8 a b and B a b 
respectively, are symmetric traceless tensors, orthogonal to u a , which we define by 

S ab = u c u d W achd (2.12) 
B ah = -\u c u d i] ac e ^ eSMl (2.13) 

where r] abcd is the covariant permutation tensor with 77 012 3 = —\J—g- 



2.1 Linearised Perturbation Equations 

Over the epoch of interest, the individual matter constituents of the universe interact 
with each other only through gravity, except for the photons and baryons (including 
electrons) whose dominant interaction with each other is via Thomson scattering of 
photons off free electrons. The variation of the gauge-invariant temperature pertur- 
bation <5t(c) along null geodesies is given by the (linearised) covariant Boltzmann 
equation |J: 

8 T (e)' + a T n e 8 T (e) = o ah e a e h + w a e a - I -9 + (1 + 4<5 T (e)) 

- o T n e (v^e a - ^^eV) , (2-14) 
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where n e is the free electron density (the effects of thermal motion of the free electrons 
is ignored), <tt is the Thomson cross section, is the baryon velocity relative to u a 
= 0(2)), and a prime denotes differentiation with respect to the parameter A 
along the null geodesic, with {u a + e a )V a A = 1. In equation (2.14) we have ignored 
the effects of polarisation. Including polarisation gives only a small correction to the 
collision term in the Boltzmann equation due to the polarisation dependence of the 
Thomson cross section. Equation (2.14) is valid for any type of perturbation (scalar, 
vector or tensor) and for any value of the spatial curvature. The evolution of the 
photon density is given by 

pin) + |0 p (7) + (3) V a g W = 0, (2.15) 

where an overdot denotes differentiation with respect to proper time along the in- 
tegral curves of u a (p^ = u a V a p^). Taking the I = 1 angular moment of the 
Boltzmann equation (2.14) gives a propagation equation for the photon heat flux: 

+ \0& + ( 3 Vtt« + |P (7 W - | (3) V.P W = a T n e (f - #)) . 

(2.16) 

Taking higher-order moments of Eq. (2.14) gives a hierarchy of equations which are 
used in the covariant numerical calculations of CMB anisotropies described in M. 

The electrons and baryons may be approximated by a tightly-coupled ideal fluid 
with energy density p( b \ pressure p^ in the rest frame of the fluid which has velocity 
u a + v£>. To linear order, the stress-energy tensor of the baryons is 

Tj? } = p^u a u b + 2(p« + p^)u {a v^ - p^h ab , (2.17) 

which shows that the baryon heat flux is (p^ + p^)v^ in the u a frame. The 
conservation of photon plus baryon stress-energy gives the propagation equations 
for the density 

and the velocity 

(P (b) + P (b Wa b) + Wa) + Hp (b) + P {b) ) dv a ] + P^f ~ ^aP^ = 

- a T n e (ip^vV - gW) . (2.19) 

In this paper we consider only scalar perturbations. In this case, the magnetic 
part of the Weyl tensor B ab and the vorticity w ab vanish identically. The electric part 
of the Weyl tensor £ ab and the shear a ab do not vanish, and satisfy the propagation 
equations 

£ ab + 0£ ab + |k[3(/o + p)a ab + 3 (3) V (a q b) - h ab ^V c q c - M ab - 6Tt ab ] = 

(2.20) 

°ab + \QOab ~ ^ ( a W b) + |/l ab (3) V C ^ C + S ab + \nil ab = 0, 

(2.21) 
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where k = 8ir, and the constraint equations 



{3) V b S ab - ±/c[2 (3) V„p + 26q a + 3^V b 7t ab ] = (2.22) 
(3) V b a ab - Pv a # - Kq a = 0. (2.23) 

The density, pressure, heat flux and anisotropic stress appearing in these equations 
are total variables obtained by summing over all matter constituents. 

For scalar perturbations, the temporal and spatial aspects of the problem may 
be separated by expanding all first-order gauge-invariant variables in tensors derived 
from the scalar harmonic functions Q^ k \ which are defined covariantly as eigen- 
functions of the generalised Helmholtz equation ( 3 )v 2 Q*- fe ) = jzpQ^ |12| satisfying 
q(k) _ 0(1). Specifically, we have 

(2.24) 
(2.25) 
(2.26) 
(2.27) 

k 

The scalar expansion coefficients, such as , are themselves first-order gauge- 
invariant variables, and they satisfy ^'V a xjf S> = 0(2). The spatial vector Q^f 1 and 
the spatial tensor Q^ b \ which is symmetric and traceless, are defined by 

Qi k) = f (3) v a Q«, q« = S (3) v (a (3) v 6) g« - \h a a k) - 

(2.28) 

Some useful properties of the scalar harmonics and derived tensors are summarised 
in the appendix to Bruni et al. This completes the definitions of quantities 

required in this paper. Further details of our notation and conventions may be 
found in 





= E^fQi fc) , 


z a 


= E k iz k Q^ ] 




k 


W ab 


k 


£-ab 


_ k 2 /K n {k) 
k 


Oab 


= E ^kQ^ab 
k 




= E4 fc) Qi fc) , 

k 


W a 


= E^i fc) - 

k 



3 A Covariant Expression for the Temperature 
Anisotropy 

The gauge-invariant CMB temperature anisotropy along a given direction is obtained 
by integrating the Boltzmann equation (2.14) along the null geodesic (whose tangent 
projects onto the given direction) through the observation point. Before integrating 
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Eq. (2.14), it is convenient to rewrite the first-order factor multiplying 1 + 45t(c) on 
the right-hand side in terms of gauge-invariant variables as follows: 

1 (7)' 1 

L 9+ P=-l— (e B (3) W W - ^V a q^) , (3.1) 
3 4pW 4pW ^ " / v 1 

where we have made use of the equation of motion of the photon density, Eq. fl2.15|) , 
and [u a + e a )( 3 )V a A = 1. At this point, we specialise to scalar perturbations and 
introduce the harmonic expansions of the gauge-invariant variables given in the 
previous section. We have 



P' 



i_(3) V a g ( 7 ) =E | g (7) QW (32) 
k 



L ea (3) v yr) = £ X M Q W (3 . 3) 



pW 



so that 



^ ( (3) v° g « - e a (3) v>«) =~E (I** - f H Q {k) - \ E(4 7) Q (fc) )'_, 



(3.4) 



where we have used the equation 



jgjM = -11^ + 11^-^', (3.5) 

which follows from taking the spatial gradient of Eq. ( |2.15| ) and harmonically ex- 
panding the result, to eliminate in favour of Z^ and the acceleration. Integrating 
the Boltzmann equation (2.14) along the null geodesic connecting the reception point 
R (where A = Xr) and a point in the distant past (where A = Aj), we find 



(Me)) R = -jE(^ (7 V fc) 
^ k 

+ E t R ^ T [i°ke a e b QW _ 1 (k Zk _ l 9w \ Q (k) + Wke a Q{ , 



/A, 



rfA 



+ E / R ~r'e- T [&7#Wg£> - fe-flW + f* fe (7) Q«] dA, (3.6) 

where (M)r denotes the value of the quantity M evaluated at the point R, and r(A) 
is the optical depth along the line of sight, defined by 

t(A) = / n e aT<i\. (3.7) 
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On angular scales larger than 8' we may approximate the visibility function —r'e~ T 
by a delta function whose support defines the last scattering surface (the instanta- 
neous recombination approximation). With this approximation, Eq. (3.6) integrates 
to 



(Me))* = E (\4 l] Q (k) + ^ } e a e b Q% - v^QS 



+ E t R {i h {usQ {k) y + iQ (k) ) - \Z k Q {k) 



f(iu*Q< fc)/ + Hw k Q^)} dX, (3.8) 



where the point A (where A = Xa) is the point of intersection of the null geodesic 
with the last scattering surface, and we have used the result 



e a e b Q {k) 



s {SQ (k)y + i Q 



(k) 



(3.9) 



We have dropped a direction independent (monopole) term evaluated at R from 
Eq. (3.8) since it will eventually be cancelled by other monopole terms in the inte- 
gral. Note that the integrand in Eq. (3.8) contains only kinematic gauge-invariant 
variables (the shear, the spatial gradient of the expansion 9 and the acceleration), 
which simplifies the next stage in the integration. In [^] we gave a more general 
expression for the anisotropy, valid for all perturbation types, but the integrand 
involved the spatial gradient of the baryon density which could only be replaced 
by the spatial gradient of the total density if the universe is baryon dominated at 
recombination. The expression (3.8) proves to be more convenient for the discussion 
of CMB anisotropics in multicomponent universes where only scalar perturbations 
are present. 

Integrating the last term in Eq. (3.8) by parts twice, we find that 



Ok (jkiSQ^J + |Q( fc )) - \Z k Q^] + f [w k QW + Hw k Q 



>(*) 



dX 



+ 



A,i 



R 
A 



k W k 



Q {k) dX. (3.10) 



The integrand on the right-hand side of Eq. (3.10) may be simplified by using the 
linearised propagation and constraint equations for the shear and the electric part 
of the Weyl tensor. The harmonic expansions of equations (2.20) and (2.21) give 



(I) 2 + + 55«P(7^fc + Qk) + |«P0(37 - 1W - 5«p7r* = 

j(<7 fe + \6o k - w k ) + (|) 2 $ fc + \npn k = 0, 



(3.11) 



(3.12) 
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where 7 is defined by p = (7 — l)p, and the constraint equation (2.23) gives 



(|) 2 [Z k - (l - ff ) Ok] + «:pg fc = 0. (3.13) 



In these equations, the scalar variables and tt^ are the harmonic expansion coeffi- 
cients of the total heat flux and anisotropic stress. They are related to the component 
variables q£' and 7r£* (defined by Eq. (2.25)) by 

pq k = J2P il) i\ pn^Y.P^^K (3-14) 

i i 

where the sums are over individual components i. 

Evaluating the integrand in Eq. (3.10) by differentiating the shear propagation 
equation (3.12), substituting for qu and Zk from equations (3.11) and ( |3.13|) , and 
using the zero-order Friedmann equation 



we find the result 



H 2 + § = \ K p, (3.15) 



(fO' + §!fa - z *) - H = -2^, (3.16) 



which is true for scalar perturbations, independent of the matter description and spa- 
tial curvature. With this, we obtain our final result for the temperature anisotropy 
(which is exact in linear theory on angular scales where instantaneous recombination 
is valid): 



(Me))« = E ([i*i 7) + f fa - *>k)] Q {k) ) A - E ([4 b) + ° k VQ { * 

k k 

E (4 7) ^Qi?) a - 2 E £ R ®kQ {k) dX, (3.17) 



T 16 



k k 



where we have dropped a (frame-dependent) dipole term evaluated at R since such 
a term cannot be distinguished from a first-order peculiar velocity of the observer at 
R. The final term in Eq. (3.17) describes the Rees-Sciama effect, which only makes a 
small contribution to the anisotropy in K = models that are matter dominated at 
recombination (for K = 0, $*: is approximately constant while a mode is outside the 
horizon and during the matter dominated era on all scales). The third term on the 
right-hand side of Eq. (3.17) represents a small contribution to the anisotropy from 
photon anisotropic stress at last scattering. The sum of the first and second terms 
dominates the CMB anisotropy in a K = universe, with the relative importance 
of each term being dependent on and H Q . Expression (3.17) is a generalisation of 
the result given by Dunsby in Section 5 of which was valid only for universes that 
are fully baryon dominated at last scattering and are spatially flat. We shall see in 
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Section [| how the result in |9j (actually a corrected version of it) may be obtained 
from Eq. (3.17) in the limit of baryon domination at recombination. A similar result 
to Eq. (3.17) is derived in [111] in terms of Bardeen's gauge- invariant variables. 

In deriving Eq. (3.17), we have not made an explicit choice for the velocity u a . 
Each of the four terms on the right-hand side is frame-independent, which follows 
from the fact that under a change of frame u a i— > u a + v a , where v a is a first-order 
relative velocity (u a v a = 0(2)), £ a b, E a b and -K ab are invariant, while v a b > h- > ti^ — v a 
and 

(Tab o ah + i3) V (a v b) - \h ab ^V c v c . (3.18) 

The frame-invariance of the right-hand side of Eq. (3.17) is necessary since St{c) 
is invariant in linear theory, up to the dipole terms that we have dropped from 
Eq. (3.17). 

The non-integral terms on the right-hand side of Eq. (3.17) are evaluated at the 
point A on the last scattering surface, which lies on a null geodesic through the 
observation point R. However, it is only necessary to locate A to zero-order since 
the displacement from the "true" position is first-order, which leads to only a second- 
order error when evaluating a first-order variable. We shall return to this point in 
Section [| where we discuss some of the gauge issues associated with the placement 
of the last-scattering surface in the standard calculations of the Sachs- Wolfe effect. 



4 CMB Anisotropy in a Universe Dominated by 
Baryons at Recombination 

In cosmological models that are baryon dominated at recombination, the covariant 
result for the temperature anisotropy, Eq. (3.17), may be cast in a more familiar 
form, which aids direct comparison with other such results in the literature (for 
example @) and physical interpretation. 

Using the propagation equations (3.12) and (3.11) for the shear and the electric 
part of the Weyl tensor, and the harmonic expansion of the constraint (2.22): 

2 (I)" i 1 i*p (X h + [l - f 1 ] vr fc ) - K P 6q k = 0, (4.1) 

we find in the limit 

p->p(*>, p-,0, X k ^4 b \ q k ^v ( k b \ n k ^0, (4.2) 
at last scattering that 



(4.3) 
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where we have added and subtracted X% /3 making use of Eq. ( |4.1| ), and 

+ + (4-4) 

It follows that in a universe which is baryon dominated at last scattering, the tem- 
perature anisotropy from scalar perturbations in the instantaneous recombination 
approximation becomes 



fc 

+ E [** + e a Qi fc) ) , - 2 E r ^Q {k) dX. (4.5) 



fc k 



The first set of terms in square brackets on the right-hand side of Eq. (4.5) give the 
"monopole" contribution (at last scattering) to the temperature anisotropy. The 
term (X^/A - X® /3)Q^ arises from entropy perturbations, which may be char- 
acterised covariantly by a vector where 

c(t&) = JL(3) V Jn) _ J_(3) V /) U 6) 

° a ~ 4pW ap p (b) VaP ' 1 DJ 

For adiabatic initial conditions, entropy perturbations vanish at last scattering on 
large scales, where tight-coupling between the baryons and photons still holds. The 
second "monopole" term, —Q k Q^/3, is the usual Sachs-Wolfe contribution to the 
temperature anisotropy [[[]], whose effect is modified on small scales by the third 
and fourth "monopole" terms. As noted by Ellis and Dunsby M, the third term is 
rarely seen in analytic calculations of the Sachs- Wolfe effect, although it is present 
in Panek's result |TD[. The omission arises from subtle gauge effects at the last 
scattering surface which we discuss in Section |5]. The final monopole term is a small 
correction arising from the non-stationarity of the potential The terms under 
the second summation in Eq. (4.5) make a "dipole" contribution to the temperature 
anisotropy, and are important on small angular scales. 

In a K = universe, Eq. (4.5) may be written in the form 



(Me)) R = E - \4 b) - & + fir 1 ** + ln^ k ] qW 
k 

+ E (fWfc 1 [** + e a Q ( a k) ) A ~ 2 E f R ^Q (k) d ^ (4-7) 



where 7i k = SH/k is the ratio of proper wavelength to the Hubble radius. This 
result corrects that given by Dunsby in Section 5 of (his equation (61); note also 
the difference in metric signature from that adopted here). Note, in particular, that 
the dominant contribution to the CMB anisotropy from adiabatic perturbations on 
large scales is — $fc/3, which is a factor of 3 smaller than the result in ||. 
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Figure 1: Contributions to the CMB temperature anisotropy from the conditions on 
the last scattering surface (z ~ 1050) in a K = universe with Qj, = 1, and Hq = 
50kms _1 Mpc -1 . Only adiabatic scalar perturbations are considered, and <&k is indepen- 
dent of k initially (so that we are plotting transfer functions). The curves should be mul- 
tiplied by the initial <]?/% to give the actual conditions on the last scattering surface. The 
solid lines are calculated from Eq. (3.17); the "monopole" contribution (negative on large 
scales) is X^ /A + S(&k — Wk)/k, and the "dipole" contribution (positive on large scales) is 
— CTfc— Vu ■ The dashed lines are given by the approximate expression (4.7); the "monopole" 

contribution (negative on large scales) is X^ /4- xf^ /3- <S> k /3 + 2H- 1 $ k /3 + 2H k 2 <S> k /9, 
and the "dipole" contribution is 27i k ~ l {<& k + H~ 1 & k )/3. The small discrepancy between 
the solid and dashed curves is due to the small but non-zero values of p*- 7 ** / p and j p at 
last scattering. 



In Fig. [I] we plot the first two summands in equations (3.17) and (4.7) as a 
function of k, with S = 1 at the present, in a K = universe, in the limit Qi> = 1, 
where Qb is the present-day baryon density in units of the critical density. We take 
Hq = 50kms _1 Mpc -1 and consider adiabatic perturbations with only the fastest 
growing mode present. At early times, we take to be independent of k. The actual 
conditions on the last scattering surface are obtained from the transfer functions of 
Fig. [I] by multiplying by the initial values of $fc (which are Gaussian random variables 
in most inflationary theories). The gauge-invariant variables on the last scattering 
surface are obtained from an accurate Boltzmann code employing covariant, gauge- 
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invariant variables, with adiabatic initial conditions |7| . The agreement between the 
approximate expression (4.7) and the "exact" expression (3.17) is good over the full 
range of k depicted. The small discrepancy between the sets of curves is due to the 
small but non-zero values of p^ / p and p^ / p on the last scattering surface (which 
is located at z ~ 1050 in this model). 

On large angular scales (small k), the "monopole" terms dominate the CMB 
anisotropy, giving the familiar Sachs- Wolfe plateau for a scale-invariant spectrum 
of initial conditions. On smaller angular scales the "monopole" term oscillates in 
k due to the coherent acoustic oscillations in the photon-baryon plasma that occur 
for modes inside the (sound) horizon. These oscillations, along with the oscillations 
in the "dipole" on small scales, determine the structure of the Doppler peaks in the 
CMB power spectrum. For adiabatic perturbations in a K = universe, the first 
zero in the "monopole" contribution to the anisotropy (the initially lower curves 
in Fig. [I]) occurs where Tik ~ V(2/3). This follows from Eq. (4.7), the fact that 
initially adiabatic perturbations are still adiabatic at last scattering on such scales 
(see Fig. |2[), and the approximate stationarity of the potential in the matter 
dominated era of a K = universe. This effect was noted recently by Ellis and 
Dunsby ||, and was also discussed by Hu and Sugiyama JUj. Ellis and Dunsby 
suggested that this zero in the "monopole" contribution should be observable as a 
zero in the CMB power spectrum on angular scales ~ 50'. However, as noted by Hu 
and Sugiyama and evident in Fig. [1], the "dipole" contribution is already significant 
on these angular scales, with its first maximum occurring close to the zero in the 
"monopole" . This tends to wash out the zero (actually a dip once the k modes are 
summed over) in the power spectrum, as shown by Hu and Sugiyama [|14]] in their 
Fig. 4. 

In Fig. ^| we plot the individual contributions to the CMB anisotropy in the 
baryon dominated limit from equation (4.7). On large scales (small k), the domi- 
nant contribution is from the usual Sachs- Wolfe term — $fc/3. The effect of entropy 
perturbations is negligible on large scales since the baryon-photon fluid is still tightly- 
coupled at last scattering on these scales, and our initial conditions are adiabatic. 
However, on small scales the entropy perturbations add coherently with the term 
2'H~ k 2 &k/9 to reduce the "monopole" contribution to the anisotropy. These terms 
along with the "dipole" term determine the structure of the Doppler peaks. 

We introduced xj® into Eq. (|4.3|) to show explicitly the contribution from entropy 



perturbations at last scattering, and to facilitate comparison with other results in 
the literature. However, this decomposition into an entropy perturbation term and 
terms involving the potential is rather unnatural, and can be replaced by an 
expression which allows a more physical interpretation of the temperature anisotropy 
on intermediate and small scales. We use Eq. ( f4.1| ) to eliminate X^ and choose u a to 
coincide with the baryon velocity [v^> = 0; if the universe is not baryon dominated 
at recombination, then the energy-frame should be employed instead of the baryon- 
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Figure 2: Approximate contributions to the CMB temperature anisotropy from the con- 
ditions on the last scattering surface (z ~ 1050) in a K = universe with flj, = 1, and 
Hq = 50kms _1 Mpc -1 . Only adiabatic scalar perturbations are considered. On large scales 
the usual Sachs-Wolfe term (dotted line) — dominates. For adiabatic initial condi- 
tions, the effect of entropy perturbations (dashed-dotted line) is only important on small 
scales where the tight-coupling approximation ceases to hold and adiabaticity is broken. 
The coherent sum of the entropy perturbations and the term 2H~^ 2 <&k/Q (solid line), as 
well as the "dipole" term 2H^ 1 (<&k + -ff _1< &fc)/3 (dashed line) determine the structure of 
the Doppler peaks. 

frame) : 

\X^ + f(** - «*) - R (7) - i*k - l§* k + g**, (4.8) 

where X^ is the harmonic expansion coefficient of the spatial gradient of the photon 
energy density in the energy-frame. This result actually holds in the energy-frame 
of any model that is matter dominated at last scattering, and so includes standard 
CDM and variants. We see that in a K = universe, which is matter dominated 
at recombination, the "monopole" contribution to the anisotropy consists of the 
Sachs- Wolfe term, —<S>kQ( k )/3, which is dominant on large scales, a term X^Q^/A 
arising from inhomogeneity of the radiation temperature in the energy-frame, which 
is important on intermediate and small scales, and the small term 2H~ 1 <&kQ^ /3 
arising from non-stationarity of the potential at last scattering. Since the last scat- 
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tering surface is well approximated by a surface of constant radiation temperature, 
the spatial gradient of the radiation energy density, in the energy frame, across the 
last scattering surface describes a distortion of the last scattering surface relative to 
the surfaces of simultaneity of the energy- frame (which are defined by the average 
motion of all matter in the universe). The distortion of the last scattering surface 
causes photons to incur extra redshift, due to the expansion of the universe, during 
propagation from last scattering to the point of observation. 

Finally, if we consider initially adiabatic perturbations in a K = baryon domi- 
nated universe on large enough scales that entropy perturbations may be neglected, 
we may replace X^ /4 by the spatial gradient of the baryon energy density in the 
baryon- frame, xf®/3. Replacing with 37i\X^ /2 (from the constraint fli.lD with 
K = 0) and neglecting terms involving we obtain the covariant equivalent of 
Panek's result [|10|j : 

(M^)* = E (f <*f [i - §w2] Q {k) ) A + E {xPhhSqw) . (4.9) 



5 Comparison with the Sachs- Wolfe Result 

It is instructive to compare the preceding analysis with the original calculation of 
Sachs and Wolfe (lj . Assuming that the radiation is nearly isotropic in the frame of 
the baryons at last scattering, the temperature of the radiation at R (in the baryon- 
frame) is given in terms of the radiation temperature at the point A in the last 
scattering surface as 



R 



(5.1) 



where the redshift z along the null geodesic connecting A and R is given by JT| 



l + z 



Sr 

Sa 



1 

2 Jr/ A 



VR 



(d v h i je l e j + 2d ri h i oe l ) dr] 



(5.2) 



where Sr = S(r] R ) is the background scale factor at R and similarly for S A , h^ v is 
the metric perturbation in the comoving gauge (with ft, o — 0), rj is conformal time, 
and e l (i = 1, 2, 3) is the direction of photon propagation in the background model. 
Differencing the temperature between two directions on the sky we find 



/AT 



V T 



R 



(f).+(¥ 



-A 



>IR 



'IE 



drjhije 1 ^ +2d ri h i0 e l ) dr\ 



(5.3) 



where (AM) E = (M) A -(M) B . The (AS) E term on the right-hand side of Eq. (|J) is 
gauge- dependent (first-order gauge transformations of the form rj i— > 1] + f(x l ) / S(rj), 
which preserve the gauge- conditions implicit in Eq. (|5.2|) , move the field S around in 
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the real universe), while the term (AT)e is gauge-invariant. It follows that the sum 
of the first two terms on the right-hand side of Eq. ( |5.3|) is also gauge- dependent. (A 
compensating gauge-dependence in the integral term preserves the necessary gauge- 
invariance of (AT)#.) This sum may be expressed in terms of the gauge- dependent 
photon density perturbation defined by 

1 f \ T- 

where = T^(rj) is the background radiation temperature. Using the constancy 
of T^S, we find 



) = \a(5^) e . (5.5) 



The significance of this result is that it is a first-order quantity. To evaluate A(S^)e 
to first-order, we need only locate the last scattering surface to zero-order, since lo- 
cating the last scattering surface correctly amounts to a first-order displacement in 
a first-order quantity. The implicit choice made by many authors is to evaluate 
Eq. (|5.5|) on the background last scattering surface (over which rj is constant). With 
this choice AS and AT' ) are both equal to zero, and trivially we obtain the result 
A(5^)e evaluated on the background last scattering surface. Note that Sachs and 
Wolfe J[| explicitly made this choice, although they also made the (gauge-dependent) 
assumption that the radiation temperature was constant on the background last scat- 
tering surface. Although it is certainly necessary to locate last scattering correctly 
to calculate A(T) E and A(S)e separately, this is not true for the difference A(TS)e, 
which is all that is required for CMB calculations. Physically, this effect arises from 
the compensation effect: the extra redshift along a given direction due to the differ- 
ence between the scale factor on the true last scattering surface and at the point of 
intersection of the geodesic with the background (zero-order) last scattering surface 
cancels out the difference in temperature between the same points on the real and 
true last scattering surface. Note that such problems do not arise in the kinetic the- 
ory calculation of Section |3], since the optical depth (which determines the position 
of the last scattering surface) multiplies only first-order variables, and so is needed 
only to zero-order. 

The term A(5^)e/A is often omitted from Sachs- Wolfe type analyses (for exam- 
ple, HT5I), effectively being absorbed into a (direction-dependent) "expected" tem- 



perature, T = TaSa/Sr (see |2|, [RJ for alternative choices). Along with the observed 

(5.6) 



temperature T R , this defines a first-order temperature variation ST/T by 

5T _T R -f 



T T 

This should not be confused with the gauge-invariant temperature fluctuation from 
the mean <5r(e), used in the rest of this paper. The quantity ST/T is gauge- dependent 
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through the gauge-arbitrariness in the scale factor S. Differencing between two 
directions on the sky, one finds 



A {ST/T) R = A / (d v h ij e l e 3 + 2d ri h i0 e l ) drj 



(5.7) 




Of course, this result is entirely equivalent to Eq. ( |5.3|) , one needs only the result 

/ATA A(TS)e , . 

\~>E~lfsj7' (5 ' 8) 

which follows from the definition of T. This appears to be the step where confusion 
sometimes arises in Sachs-Wolfe type analyses, with authors comparing the gauge- 
dependent 5T/T with the results of observation, instead of the physically relevant 
(AT/T)r (or 5t{g) which is easily derived from the former). As noted in 0, it is 
not enough just to difference 5T/T; in general, the result will not be gauge-invariant 
and will not include the contribution from 5^ at last scattering. 

Finally, we demonstrate how the gauge-invariant contribution to the 

anisotropy arises in the Sachs- Wolfe calculation. If we approximate the (vorticity- 
free) baryon motion as being geodesic through recombination to the present, we 
can choose a comoving gauge with hio = (the synchronous gauge, in which the 
surfaces of constant r] are surfaces of simultaneity for the baryons). The result given 
in Eq. ( |5.3| ) (with hio = 0) still holds in this gauge, and the conformal time r\ satisfies 
V aV = u^>/S(t]), where u$ is the baryon velocity. It follows that the background 
scale factor satisfies ^V a S = and S = d v S/S = 9S/3 + 0(1), where 9 = V a u^ 
is the covariant expansion of the baryon-frame. With this gauge-condition, we have 
removed the gauge-freedom to perform the transformation r] i— ► rj + f(x l )/S (the 
gauge- condition ^V a r/ = forces f(x l ) to be a constant), with the result that AS/S 
is gauge-invariant. Under these conditions, we may express the first two terms on 
the right-hand side of Eq. ( |5.3j ) in terms of covariant, gauge- invariant variables as 
follows: 

1 f A ^" M + ieu?)^ 



AJb \ pW 3 
= EM|^Q( fc) )^ M) 

k 

where A and B are the points of intersection of the two geodesies through the 
observation point R with the last scattering surface, dx a lies in the last scattering 
surface (dx a u^ = 0(1)), and we have replaced 3S/S by the covariantly defined 9 
which is correct to the required order. In this manner, we recover the 4 (7 V fc )/4 
contribution to the temperature anisotropy. 
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6 Conclusion 



Starting from a covariant and gauge-invariant formulation of the Boltzmann equa- 
tion, we have derived a new expression for the CMB temperature anisotropy under 
the instantaneous recombination approximation, valid for scalar perturbations in 
open, closed and flat universes. Our expression uses only covariantly-defmed vari- 
ables, and is manifestly gauge-invariant. The result is more useful in multicomponent 
models with scalar perturbations than earlier covariant results || 0. In the case 
of a universe which is baryon-dominated at recombination, we find a simple expres- 
sion for the anisotropy which corrects a similar result by Dunsby ||. By making 
use of numerical solutions to the perturbation equations, we have discussed the con- 
ditions on the last scattering surface and their contributions to the characteristic 
features of the CMB power spectrum. We ended with a discussion of the original 
Sachs- Wolfe calculation for the temperature anisotropy. We have discussed why it 
is not necessary to locate accurately the last scattering surface in such calculations 
(because of the compensation effect), and how the extra term in the Sachs- Wolfe 
calculation, reported recently by Ellis and Dunsby ||, is missed in many calcula- 
tions which employ a gauge-dependent "expected temperature", since the angular 
dependence of this temperature is often overlooked. For a universe which is matter 
dominated at recombination, but not necessarily adiabatic, the extra term is the 
spatial gradient of the radiation energy density in the energy- frame, (X^Q^ /4) A . 
For models with isothermal surfaces of last scattering, this inhomogeneity describes 
a distortion of the last scattering surface relative to the surfaces of simultaneity 
of the energy-frame. The extra redshift incurred by this distortion is a significant 
component of the temperature anisotropy on intermediate and small scales. 
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